function OpenCRTBP_SE_BaseFigure
GM_sun=1.327124400*10^(11);     %Sun
GM_earth=3.986004418*10^(5);    %Earth
GM_moon = 4.9048695*10^(3);     %Moon
%[km^3/s^2] Gravitational Parameters
r12_se = 1.4959965e+8;      % Distance between sun and the earth
TU_se = 5.0226757e+6;       % Time unit
LU_se = r12_se;             % Length unit
VU_se = LU_se/TU_se;        % Velocity unit
u_se = GM_earth / (GM_earth+GM_sun);
u_em = GM_moon / (GM_moon+GM_earth);
%[] Mu constant
%    se --> Sun-Earth system
%    em --> Earth-Moon system
re = 6378/r12_se;
rs = 696340/r12_se;
figure
hold on
[L1,L2,L3,L4,L5] = librationPoints(u_se);
Lse = [L1,L2,L3,L4,L5];
for ii = 1:5
    plot3(Lse(1,ii),Lse(2,ii),Lse(3,ii),'bo','MarkerFaceColor','b','MarkerSize',2)
end
[xs,ys,zs] = ellipsoid(-u_se,0,0,rs,rs,rs,50);
[xe,ye,ze] = ellipsoid(1-u_se,0,0,re,re,re,50);
Earth = imread('Earth.bmp');
surface(xe,ye,-ze, ...
        'CData',Earth, ...
        'EdgeColor','None', ...
        'FaceColor','TextureMap')
surface(xs,ys,zs,'FaceColor',[255 165 0]/255,'EdgeColor','None');
xlim([-1.3,1.3]);
ylim([-1.3,1.3]);
axis square
set(gca,'FontSize',20) 
view(3)
end